arXiv: 1502.06135v2 [hep-th] 3 Dec 2016 


Preheating with Fractional Powers 

Hossein Bazrafshan Moghaddam, Robert Brandenberger 
Department of Physics, McGill University, Montreal, QC H3A 2T8 Canada 

We consider preheating in models in which the potential for the inflaton is given by a fractional 
power, as is the case in axion monodromy inflation. We assume a standard coupling between the 
inflaton field and a scalar matter field. We find that in spite of the fact that the oscillation of the 
inflaton about the field value which minimizes the potential is anharmonic, there is nevertheless 
a parametric resonance instability, and we determine the Floquet exponent which describes this 
instability as a function of the parameters of the inflaton potential. 


INTRODUCTION 

There has been a lot of recent interest in large field models of inflation in which the potential is given by a fractional 
power of the field. An example which is currently attracting much attention is axion monodromy inflation (see [II 
for the initial paper and [2j for a recent review). Axion monodromy models are attractive since they may provide a 
natural realization of large field inflation in the context of superstring theory. Large field models of inflation EIH 
are advantageous since in such a context the slow-roll trajectory is a local attractor in initial condition space [5], even 
including metric fluctuations [B]. In contrast, for small field models the initial velocity of the inflaton field needs to 
be fine tuned, thus creating a potential initial condition problem [7]. As is well known [8], from the point of view of 
observations, large field models of inflation are interesting since they may lead to a significant tensor to scalar ratio. 

There are challenges to obtain large field inflation models. For field values |</>| > m p i corresponding to large field 
inflation ( m p i is the Planck mass), there is the danger that gravitational corrections will lift the potential and prevent 
slow rolling of the field, unless there is a symmetry such as shift symmetry (5] which protects the small mass required 
for large field inflation. In string theory, there is a further challenge of obtaining large field inflation: we expect 
the field range of the candidate inflatons, e.g. the moduli fields or the fields associated with brane separations, to 
be small, and hence incompatible with large field inflation. Monodromy inflation m provides a possible resolution 
of this problem, and axion monodromy is currently regarded as the most promising implementation of the idea of 
monodromy inflation in the context of string theory [2]. For this reason, there has been a lot of recent activity on 
this topic (see e.g. mm for a selection of recent papers). 

To our knowledge there has been very little work on reheating in models in which the inflaton potential is given by 
a fractional power of the field <f> (but see [23] for a study of reheating in the context of the initial brane monodromy 
inflation model of HD]). The goal of this paper is to consider preheating in models in which the inflaton potential is 
given by a fractional power of the field, and in which the inflaton field is coupled to a second scalar matter field \ via 
a usual interaction term c f >2 x 1 2 IOur goal here is to study the changes in preheating which arise due to the fact that 
the oscillations of the inflaton about the field value which minimizes the potential are not harmonic. 

Reheating is an integral part of any inflationary model. Without reheating the universe after inflation would be 
empty and cold. Reheating was initially discussed using first order perturbation theory [27i - B5] . In [BB] and [BTj it 
was realized that an oscillating inflaton field will induce a parametric resonance instability in fields which couple to 
the inflaton (and also in the equations for the inflaton fluctuations themselves). This instability was then worked out 
in detail in [32HM] and given the name “preheating” 2 . 

Reheating has been analyzed mostly for inflaton potentials of the form V(4>) = ^m 2 </> 2 and V(<f>) = A <f> 4 , and for a 
matter field \ (treated for simplicity as a scalar field) coupled to (f> via an interaction Lagrangian £ int ~ </> 2 \ 2 - For 
both types of inflaton potential, the linearized equation of motion for fluctuations of \ leads (neglecting the expansion 
of the universe) to a harmonic oscillator equation with a periodically varying mass. In the first case of a massive 
inflaton field, the equation for Fourier modes of \ is an equation of Mathieu type [23, in the second case, the case of 
a massless inflaton field with quartic potential, to a generalized equation. In both cases, the equation of motion for 
the Fourier modes \k is of Floquet type, and the general theory [36] tells us that the modes will evolve as 

Xk ~ Fi, fc (t)e /ife ‘ + P2,kit ) e ~ flkt , (1) 


1 For a perturbative study in the context of axion monodromy inflation see m- 

2 For recent reviews of reheating the reader is referred to m and m■ 
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where Pi,k(t) (i = 1,2) are periodic functions, and pk is a non-negative exponent called the “Floquet exponent”. 
Values of k for which > 0 form resonance bands. 

Past work on reheating has shown that the efficiency of the parametric resonance instability depends on which 
class of model one is considering. We are interested in inflaton potentials with fractional power of the field. As in 
the canonical inflation models, the inflaton will be oscillating about its ground state. However, the oscillation is not 
sinusoidal. Given the recent interest in axion monodromy inflation it is hence of great interest to study the efficiency 
of reheating. This is the problem we address in this paper. We will study the structure of the resonance bands 
and estimate the value of the Floquet exponent within the resonance bands. Our results show that the parametric 
instability is of broad resonance type and hence efficient. It will lead to the transfer of the inflaton energy density from 
the homogeneous inflaton background to matter fluctuations within less than a Hubble expansion time. Our result will 
allow a more precise determination of the parameters of models such as axion monodromy inflation from cosmological 
observations. As has been emphasized recently (see e.g. [371 d]), for a given action the values of the cosmological 
parameters such as the scalar spectral index n s and the tensor to scalar ratio r depend on the reheating history. One 
of the reasons for this is that changing the duration of reheating will change the current physical wavelength which 
exits the Hubble radius a fixed number of e-foldings before the end of inflation. Our result also has implications for 
baryogenesis. Since the energy transfer to regular matter after inflation is rapid, the “reheating temperature” is high 
and high-scale baryogenesis mechanisms become feasible. 

In the next section we give a brief review of axion monodromy inflation as an example of the type of models we 
have in mind. Section 3 is the main one in which we study the equation of motion for the fluctuations of a matter field 
coupled to the inflaton and establish that preheating is efficient. We work in the approximation that the expansion 
of space is negligible. Since the instability in our model is of broad resonance type and not of narrow resonance, 
including effects of the expanding space will not change the qualitative features of our results. I The most important 
change is that instead of exponential growth of the fluctuations we will have exponential growth damped by a power 
of the cosmological scale factor In Section 4 we show that this approximation is self-consistent and we analyze various 
back-reaction effects. 

We use standard notation in which a(t) is the cosmological scale factor, t is physical time, and H(t) is the Hubble 
expansion rate. 


AXION MONODROMY INFLATION 

The low energy limit of string theory provides many scalar fields which can be candidates for an inflaton field (see 
e.g. iil9| for a recent comprehensive review of inflation in the context of string theory). Examples are radii of extra 
dimensions or distances between brane-antibrane pairs. Typically, however, the range of field values is small and does 
not admit super-Planck values required for large-field inflation. 

Monodromy m is a simple solution to this problem. A string connecting brane-antibrane pairs can be wrapped 
many times about an internal cycle, thus yielding large field values even if the basic field range is small. As was 
shown in T|, suitably wrapped branes induce trans-Planckian field ranges for the associated axion fields. In string 
compactifications, axions arise from integrating gauge potentials over cycles in the internal manifold. The potential for 
the axions induced by the presence of the branes of the axion field is not periodic but grows without bound, however 
with a power which is typically less than 2 3 . One example studied in [T has a linear potential, but potentials 
with fractional powers p , e.g. p = 2/3, also arise (see e.g. [10]). The potentials have oscillatory terms (caused by 
instanton effects) and thus can give interesting observational signals in the spectrum of cosmic microwave background 
anisotropies (see e.g. [40IH3] b However, for powers p < 2 which we consider here the oscillatory terms are negligible 
for field values close to the minimum of the potential (which we take to be at </> = 0). It is the dynamics near <j> = 0 
which determines the efficiency of reheating, and hence the oscillating terms in the potential will have a negligible 
effects on the reheating efficiency. We will hence here consider simple power law potentials. 


3 Very close to the origin the inflaton potential can be quadratic EH. but that field range is negligible in terms of setting up the overall 
dynamics of <j>. As we will see, what matters for the efficiency of preheating is the amplitude of <j) near the minimum of the potential, 
and this is set by the asymptotics of the potential for field values where inflation ends. We thank E. McDonough for alerting us to this 


issue. 



3 


EVOLUTION OF FLUCTUATIONS AFTER INFLATION 

As discussed in previous section, we will consider the effective potential for inflaton field to be 

VefM |0| p , (2) 

where m is a parameter with units of mass. We will also assume a simple coupling of the form 

Ant = (3) 

between the inflaton 0 and a matter field y. Such a coupling arises at tree level (but at second order in the axion 
coupling constant) between the axion and gauge fields (which we can then represent in a simplifying way by a scalar 
matter field y). Therefore in the rest of the paper we will work with the potential 

y) = rn A ~ p |0| p + * g 2 4> 2 x 2 • (4) 

If we take the matter field x t° be massless on the scale of inflation, and if we neglect self-interactions of y, then 
each Fourier mode Xk of y will evolve independently with an equation of motion 

Xk + 3Hxk + (-^2 + g 2 Ht) 2 )xk = 0 . (5) 

a(f ) 2 

This differential equation has a time-dependent frequency and a gravitational damping term. The time dependence 
is given by evolution of both the scale factor a(t) and of the background inflaton field 0(f). 

Following E03, we will neglect the expansion of the universe for the moment, and we will later show that this is a 
self-consistent approximation during the early stages of preheating) and thus take a(t) = 1 4 . Then we can rewrite 
Eq. ([5]) as 

Xk + f{t) 2 Xk = 0 , (6) 

where 

f(t) 2 = k 2 + g 2 m 2 - (7) 

If f 2 ( t ) is a periodic function of time, then the equation is of Floquet type [35j 36] and we hence there is the possibility 
of a parametric resonance instability for the y fluctuations. The strength of the resonance depends sensitively on the 
dynamics of 0 near the minimum of the potential. We will consider the case in which the homogeneous value of the 
y field is zero. Since the effective potential for the inflaton in this case is given by ([ 2 ]), the equation of motion for the 
background inflaton field is 

<j> + 3H<fi + pm 4 ~ p |0| p_2 0 = 0. (8) 

As mentioned above, we will neglect the expansion of space and later check for the self consistency of this approxi¬ 
mation. In this case, the equation becomes 

4> + pm 4 ~ p \4>\ p ~ 2 4> = 0. (9) 

We consider values 0 < p ^ 2 for the exponent in the potential. In spite of the fact that this equation is nonlinear and 
does not describe harmonic oscillation of inflaton field, the evolution of the inflaton field is periodic (periodic motion 
about the minimum of the potential). 

As we know from previous studies of preheating [34j , the efficiency of preheating is determined by the dynamics of 
(j) near the minimum of its potential, i.e. by the value of </> at the time t c when 0 = 0. This value can be determined 
by energy conservation. The energy density at the end of inflation (t = t en d) is given by 

Pend, = ^ 2 (t end ) + m 4 ~ P \(j)\ P (t end ) (10) 


4 The expansion of space can be taken into account using the methods of |331134| . 
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from which we can find 


0(*c) \J^Pend • 

Inflation ends when the slow-roll conditions break down. These are 

i</> 2 < V{<j>) and 
3 Hj> < V', 

where the prime denotes the derivative with respect to (j>. The first condition breaks down once 


(ii) 


( 12 ) 

(13) 


(14) 


For 1 < p < 2 the second condition breaks down for a smaller value of \<j>\. For 0 < p < 1 the number p in the 


numerator of (14) is replaced by y/p (2 — p). We will denote the end point of inflation by </> e and will write 


t>e — CXTflpi 


(15) 


where the value of the constant a can be read off from (14), and we have used G = m pl . It is convenient to parametrize 
the end point of inflation in terms of a dimensionless number a given by 


am = am p i. 


(16) 


This will simplify some of the later equations. 

Now we turn to the study of x particle production, i.e. to the study of the structure of the resonance in the Eq. 
Q. While the the adiabaticity condition 


dm 

dt 


< .fit ) 2 


(17) 


is satisfied then using the WKB approximation we conclude that the state tracks the instantaneous vacuum and we 
find the solution of the from 


Xk 


a k 

7W) 


exp (—i 


fm) + VW) eMi 


f(t)dt ), 


(18) 


where ak and fdk are coefficients determined by the initial conditions. 

For small field values, the adiabaticity condition © is violated. To estimate the range of field values for which 
this is the case, we can replace <j>(t) by namely 


4> — 2 m A ~ p (aMpi )p = 2a p ^ 2 m 2 , 


(19) 


since for small field values </> is constant to a good approximation. Making use of the formula ([ 7 ]) for f(t) we see that 
adiabaticity violation arises for field values obeying 

2 g 2 a^m 2 (t> > {k 2 + g 2 (j) 2 f l ‘ 2 . (20) 

Considering modes in the far infrared (i.e. setting k = 0) we find 


H < 


\/2a* to 

~vf~ 


( 21 ) 


We denote the resulting value of \(j>\ by </> max . Thus we conclude that for small values of k the condition for violation 
of adiabaticity is 


_V2am < 

Va ~ 


\[2 a 4 m 

V9 


(22) 
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In order to estimate the energy transfer via this parametric resonance instability it is important to estimate the 
width of the instability band. To a first order, this is given by the range of k values for which the approximation 
made above of neglecting k in the adiabaticity condition violation calculation is self-consistent, namely value of k for 
which k < g4>max- This gives 


— 


r P/ 4 r 


(23) 


To put this value in context, note that for p < 4 this value is parametrically larger than H , which means that the 
modes which dominate the phase space of modes which undergo parametric resonance have wavelength smaller than 
the Hubble radius H~ x at the end of inflation, and that hence neglecting metric fluctuations is well justified (see [49] 
for a discussion showing that metric fluctuations have a negligible effect if the wavelength is smaller than the Hubble 
radius). Note that even modes with larger values of k feel the resonance, but for a shorter time interval, namely a 
time interval corresponding to 


\(t>\ < g k. 


(24) 


Now we return to Eq. © and want to find amplification of Xk modes after passing once through the range of 4> 
which allows resonance. Using Eq. (19) for the range (22) of field values where resonance can take place we can write 
(by rescaling time) 


Rewriting Eq. ([h]) in terms of new time variable 


~ 2 a 2 m~t. 


= (2g) 2 a i mt 


we get 


where the rescaled wavenumber q is 


Xk + (<7 2 + T 2 )\k = 0, 


(25) 


(26) 


(27) 


q = 


VaY- 


(28) 


Equation (27) describes scattering from a negative parabolic potential in quantum mechanics. According to the 


usual analogy where the classical variable x becomes the quantum mechanical wave function, and the time variable 


in classical dynamics becomes the spatial variable in the quantum problem, the equation (27) corresponds to the 
time-independent Schrodinger equation 


V 2 V> + {E- V)ip = 0, 


(29) 


for a wave function if) in a parabolic potential V = t 2 and with energy E. This analogy has been used in the context 
of reheating in [34], and in the case of reheating in the presence of noise in [50] El]. The time range in the reheating 
problem when the field is in the range (22) corresponds to values of x where quantum tunneling occurs. We will use 
the standard methods from quantum mechanics to solve the problem. 


Well before </> enters the range (22) the solution is of the form (18),namely 


eM ' / f(T)dT) • 


Long after it passes through that region, the solution again is of the same form, namely 


X.k 




ter pH / 


f(r)dr) 




exp(* / /(r)dr). 


(30) 


(31) 


The coefficients in Eq. (30) come from initial conditions. We can find the coefficients in Eq. (31) using amplitudes of 


reflection Rk and transmission Dk for the scattering. These coefficients are given by 

Rk = — - 


le 




Vl + e^ 2 


(32) 
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and 


Dk 


e ~i<Pk 

Vl + e ~ vq2 ’ 


(33) 


where 


\R k \ 2 + \D k \ 2 = 1 


(34) 


and tpk is given by 

ifi k = ar g r( 1 + ?g ) + ^~(l + ln4)- ( 35 ) 

2 2 q z 

We also note that the dependence on the phase of incoming wave during resonance the amplitude can increase or 
decrease. Making use of these results, and including the phase 0° of the incoming wave in the calculation, we obtain 
the following matrix relation between the coefficients of the positive and negative frequency modes before and after 
the scattering: 

( a k\ — {Vl + e~' nq 2 e llfik ie~^'> q2+2l9 k \ fa® 

[ft) ~ f)® 2 - 2 ^ Vl + e-^e -^J 

Returning to our reheating problem, we see that | /3 k | 2 gives the occupation number of particles with momentum 
k after one passage of (f> through the minimum of its potential, given that the initial occupation number is |/lg| 2 . 
Considering n\ = \fi k \ and the same relation for incoming particles we get 

= exp(—7T q 2 ) + (1 + 2exp(-7rg 2 ))?r£ - 2exp(-^g 2 ) v /l + exp(-7rg 2 ) yjn° k (l + n° k ) sin(6> t ° ot ) (37) 

where 

9tot = 20 k ~Vk + arg Pi - arg a° k . (38) 

The above relation between initial number density and final number density can be applied in sequence to all time 
intervals during which the adiabaticity condition is violated. Therefore, for the j’th scattering event we obtain (in the 
case of large occupation number) 

n{ +1 ~ (1 + 2exp(-7rg 2 ) - 2sm(0 J tot ) exp (-|< ? 2 )yjl + exp(-7rg 2 ))?4 . (39) 

Defining the “Floquet exponent” /i). for the j’th scattering through 

n fe +1 = exp(27r n J k )n J k (40) 

we obtain 

/4 = ln(l + 2exp(-7rg 2 ) - 2 sin(6>J ot ) exp(-|g 2 ) y/l + exp(-7n? 2 )). (41) 

In Figure 1 we sketch the value of /K max as a function of k. 

Summing over all scattering events, considering the initial conditions n k = 0, = l,/3g = 0, and taking the phases 

0i to be randomly distributed we obtain 

nk(t) = ^exp(27r^/4) (42) 

j 

We can rewrite the exponent as 

Vfc = VeffN 
j 



(43) 
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FIG. 1: Sketch of the Floquet exponent from Eq. (41) (vertical axis) as a function of the value of k (horizontal axis) for typical 
value of sin# = 0, p = 1, a = 0.1, m = 10 _4 m p and g = 10 -3 . The values on the horizontal axis are in Planck units. We also 
mark the value of H on this axis. 


where N is the number of oscillations of inflaton field. At a time t after the end of inflation, the number N is given 

by 


N = 


t 

T /2 


(44) 


where T is the period of infalton field oscillations is T and we devided by 
inflaton field passes through the non-adiabatic region two times. 

The period T can be estimated taking <fi to be the constant value given 


two since in each period of oscillation, the 


in (191 and solving 


Tcj) = 4am , 


(45) 


since one quarter of the integrated field range over one period is am. Hence 



m 


(46) 


Therefore 


XX ~ *• ( 4? ) 

a 2 
3 


The next step is to find the value for ii e f f . Considering small values for q and averaging over sin 9tot , we find from 
(41) 

fj-eff - T ln (3) • (48) 

To find the energy density p x of the particles produced during preheating we need to integrate Ukk over all momenta 
which experience parametric resonance, i.e. with \k\ < k max 


r ‘max ^ 

p x {t) = / d 3 k'k'-exp{2Trp , eff —- F ) 

Jo 1 v 2 

- P(ln(3)-^-) 

9 a 2 


(49) 

(50) 


where k max is given by (231. Inserting the value for k max we find 

16?r 


2 mt 


Px( l ) ~ -rT9 m aVexp(M3)^r), 

9 a 2 


(51) 












which is to be compared with the inflaton density at the end of inflation which is 


p{t e ) ~ 2 m 4 a p 


( 52 ) 


We conclude that provided that the resonance lasts so long that the exponential growth factor can overcome the factor 
i| -g 2 by which (|5l|) is suppressed compared to (52), then the parametric resonance instability will be sufficiently 
strong to drain alTof the energy stored in the inflaton field at the end of inflation. Whether this is the case or not 
depends on how soon back-reaction effects take over and how soon the approximations we have made cease to be 
self-consistent. 


BACK-REACTION EFFECTS 


In the analysis of the previous section we have neglected the expansion of space. In toy models such as those studied 
in jUJEI! this is a self-consistent approximation. We must check whether neglecting the expansion of space is justified 
in the present model. First of all, note that our analysis shows that the non-adiabatic evolution of the solution of the 
equation of motion for y takes place in narrow time intervals At when <f> is close to <f> = 0. Combining (19), (22) and 
the value for H given by the potential energy density at the end of inflation, i.e. at <j> = (f> e we obtain 


H 1 m p i 

At m ’ 


(53) 


and hence we conclude that the cosmological expansion does not affect the study of y particle production during the 
resonant time interval about any one <j> = 0 crossing point. Since the period T of oscillation obeys 


TH~\t e ) < 1, 


(54) 


neglecting the expansion of space is also justified for any one given oscillation period. Finally, taking as an example 
p = 1 and m = 10~ 4 m p i, which yields a Hubble expansion at the end of inflation close to the upper bound from 
gravitational radiation constraints, we find 


mH 1 (t e ) 


v -1 


10 . 


(55) 


Thus, over one Hubble expansion time the exponential increase in p x from ( |51[ ) is sufficient to overcome the small 
factor of g 2 in (51), taking g A ~ 10~ 12 such that the corrections to the potential for <j> due to x loops does not dominate 
over the original V((j>) at the field value </> ~ m p i. Thus, we conclude that it is a reasonable approximation to neglect 
the expansion of space. 

We will now consider some back-reaction effects which could terminate or even prevent preheating. First of all, the 
fluctuations of % which are generated in the preheating instability will lead to correction terms both in the equation 
of motion for <f> and in that of Xk- We need to determine the length of time these effects can be neglected. Secondly, 
fluctuations in the inflaton field themselves could be amplified and then back-react both in the equation of motion for 
<f>, shutting off the oscillations which drive preheating, and also in the equation of motion for y, providing corrections 
to the mass term which will prevent the violation of adiabaticity which is required to obtain the resonance. 

We first discuss the possible amplification of fluctuations of 4>. The equation of motion for the 5(p k mode is 


5<p k +3HS(p k + [^— + p(p- l)m 4 p 4> p 2 }5(t>k = 0. (56) 

cr 

To solve this equation and study possible resonance effects, we can for the moment neglect the expansion of the 
universe. Since we are interested in large scale modes we can also drop the k 2 term. Rewriting the equation and 
multiplying both sides by (p 2 ~ p gives 

(j> 2 ~ p 8(j) +p(p — l)m 4 ~ p S(/) = 0, (57) 

which (using |l5| ) and ©) can be re-written as 


A^-PScp + A 2 8(j) = 0 , 


( 58 ) 
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where 


Ai = [2{otm.pl)*m*] 2 p 

A 2 = p{p — l)m 4 ~ p . 


The solutions of (58) are Bessel functions of the first and second kind and take the following form 


6 $ 


r , T/ l 2 

ClVtJly-, ~ 
V P 



t p/2 ) + c 2 VtY(-,~ 

p p 



( 59 ) 


(60) 


The solutions for S(j> have a mild time dependence (much weaker than exponential). Hence, there is no parametric 
resonance of 4> fluctuations, and thus no important back-reaction effects of these modes. 

Thus, we now turn to the discussion of the back-reaction effects of the y fluctuations. These fluctuations will lead 
to an induced potential for the inflaton in Eq. |8]), the magnitude of which is given by g 2 x 2 4 >2 ■ This term will have a 
major effect on inflaton dynamics when the interaction term become as important as the bare potential term m A ~ p (j) p . 
Therefore the criterion for this back-reaction effect to be negligible is 


Xeff « 


n i-p 


■kP-2 


(61) 


where the value of the inflaton field to be inserted is the one at the end of inflation. 

Due to nonlinearities in the x sector (which will inevitably be present due to renormalization considerations) 
there is a term in the potential of the form | Ay 4 which will add a term to the mode equation |5j) . In the Hartree 
approximation, this term will be an extra mass of magnitude A < y 2 >. Computing this effective mass by averaging 
over the modes which undergo parametric amplification, and comparing with the amplitude of the mass term which 
is driving the parametric resonance, we find that the induced mass term is negligible if 


Xeff < ^- 



(62) 


where once again the value of <t> at the end of inflation is to be used. 

The observed amplitude of density fluctuations constrains the value of the coupling constant g 2 : at one loop level 
the interaction between y and </> induces a quartic self coupling term of (f>. The coefficient of this term is proportional 
to g 4 but it is constrained to be smaller than 10“ 12 . Hence the constraint on g 2 is g 2 < 10~ 6 . In turn, the same loop 
effects induce the self coupling of y, and we expect A ~ g A . Since the value of m is also constrained to be small from 
the observed amplitude of the power spectrum of density fluctuations, we expect the first criterion (61) to be more 
stringent than the second one (62). 


To see what the criterion (611 implies, note that when Xeff reaches the limit given by (61), then the energy density 
in the g 2 <p 2 y 2 term in the Lagrangian is of comparable magnitude to the initial energy density in the inflaton field, 
and (for values of p < 2) the energy density from the nonlinear terms in the potential for y is larger. Thus, we 
conclude that the back-reaction effects studied here do not terminate the preheating instability before a fraction of 
order 1 of the initial inflaton energy density has been drained from the coherent oscillation of <j>. Hence, we conclude 
that the parametric resonance instability is very effective. 


CONCLUSIONS 

In this paper we have studied the effects on preheating which result when considering the “non-standard” inflaton 
potentials in which the potential scales as a fractional power of the field, as is the case in axion monodromy inflation. 
We have shown that the parametric resonance instability is expected to be very effective 5 . This result will allow a 
more precise determination of cosmological parameters in a given inflationary model of the type we consider, and it 
also implies that high scale baryogenesis scenarios are possible. 


5 While our work was being written up, a paper m appeared which also studied the dependence of the efficiency of preheating on the 
shape of the inflaton potential and concluded that changes from the canonical shape do not decrease the efficiency of the loss of energy 
from the inflaton condensate. 
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In this work we have used a standard 4> 2 x 2 coupling between the inflaton and the preheat fields. Whereas such 
a coupling will arise at higher orders in the coupling constant in actual axion monodromy inflation models, in such 
models there are other couplings which are more important. In work in progress we are studying preheating for these 
more realistic couplings. In these more realistic models we can then also study the possible preheating of entropy 
modes of the metric fluctuations. Such metric instabilities could lead to important constraints on axion monodromy 
inflation models [55U55] , 

Our work has been based on analytical approximations and we have not made use of the various numerical codes 
which have been developed |57ll63j . Since we are interested in the onset and overall efficiency of the parametric 
resonance instability, an analytical approach is not only sufficient but maybe superior since it provides analytical 
insight into why we obtain the dynamics which we see. 


ACKNOWLEDGEMENTS 

We thank A. Abolhasani, Y. Cai, E. Ferreira, H. Firouzjahi, E. McDonough, G. Moore, O. Ozsoy, Y. Wang and S. 
Watson for useful discussions. This work was supported in part by a NSERC Discovery grant and by funds from the 
Canada Research Chair program (RB). 


[1] L. McAllister, E. Silverstein and A. Westphal, “Gravity Waves and Linear Inflation from Axion Monodromy,” Phys. Rev. 
D 82, 046003 (2010) arXiv:0808.0706 [hep-th]]. 

[2] E. Silverstein, “Les Houches lectures on inflationary observables and string theory,” arXiv: 1311.2312 [hep-th]. 

[3] A. A. Starobinsky, “A New Type of Isotropic Cosmological Models Without Singularity,” Phys. Lett. B 91, 99 (1980). 

[4] A. D. Linde, “Chaotic Inflation,” Phys. Lett. B 129, 177 (1983). 

[5] R. H. Brandenberger and J. H. Kung, “Chaotic Inflation as an Attractor in Initial Condition Space,” Phys. Rev. D 42, 
1008 (1990). 

[6] H. A. Feldman and R. H. Brandenberger, “Chaotic Inflation With Metric and Matter Perturbations,” Phys. Lett. B 227, 
359 (1989). 

[7] D. S. Goldwirth and T. Piran, “Initial conditions for inflation,” Phys. Rept. 214, 223 (1992). 

[8] D. H. Lyth, “What would we learn by detecting a gravitational wave signal in the cosmic microwave background 
anisotropy?,” Phys. Rev. Lett. 78, 1861 (1997) hep-ph/9606387 . 

[9] M. Yamaguchi and J. Yokoyama, “New inflation in supergravity with a chaotic initial condition,” Phys. Rev. D 63, 043506 
(2001) hep-ph/0007021 . 

[10] E. Silverstein and A. Westphal, “Monodromy in the CMB: Gravity Waves and String Inflation,” Phys. Rev. D 78, 106003 
(2008) arXiv:0803.3085 [hep-th]]. 

[11] T. C. Bachlechner, C. Long and L. McAllister, “Planckian Axions in String Theory,” arXiv: 1412.1093 [hep-th]; 

[12] T. Higaki and F. Takahashi, “Axion Landscape and Natural Inflation,” arXiv:1409.8409 [hep-ph]; 

[13] R. Blumenhagen, D. Herschmann and E. Plauschinn, “The Challenge of Realizing F-term Axion Monodromy Inflation in 
String Theory,” arXiv:1409.7075 [hep-th]; 

[14] I. Ben-Dayan, F. G. Pedro and A. Westphal, “Towards Natural Inflation in String Theory,” arXiv: 1407.2562 [hep-th]; 

[15] S. Franco, D. Galloni, A. Retolaza and A. Uranga, “Axion Monodromy Inflation on Warped Throats,” arXiv: 1405.7044 
[hep-th]; 

[16] L. McAllister, E. Silverstein, A. Westphal and T. Wrase, “The Powers of Monodromy,” JHEP 1409, 123 (2014) 
arXiv: 1405.3652 [hep-th]]; 

[17] M. Arends, A. Hebecker, K. Heimpel, S. C. Kraus, D. Lust, C. Mayrhofer, C. Schick and T. Weigand, “D7-Brane Moduli 
Space in Axion Monodromy and Fluxbrane Inflation,” Fortsch. Phys. 62, 647 (2014) arXiv: 1405.0283 [hep-th]]; 

[18] R. Kappl, S. Krippendorf and H. P. Nilles, “Aligned Natural Inflation: Monodromies of two Axions,” Phys. Lett. B 737, 
124 (2014) |arXiv: 1404.7127 [hep-th]]; 

[19] F. Marchesano, G. Shiu and A. M. Uranga, “F-term Axion Monodromy Inflation,” JHEP 1409, 184 (2014) arXiv: 1404.3040 
[hep-th]]; 

[20] K. Harigaya, M. Ibe, K. Schmitz and T. T. Yanagida, “Chaotic Inflation with a Fractional Power-Law Potential in Strongly 
Coupled Gauge Theories,” Phys. Lett. B 720, 125 (2013) arXiv:1211.6241 [hep-ph]]; 

[21] J. P. Conlon, “Brane-Antibrane Backreaction in Axion Monodromy Inflation,” JCAP 1201, 033 (2012) arXiv: 1110.6454 
[hep-th]]; 

[22] X. Dong, B. Horn, E. Silverstein and A. Westphal, “Simple exercises to flatten your potential,” Phys. Rev. D 84, 026011 
(2011) arXiv:1011.4521 [hep-th]]. 

[23] R. H. Brandenberger, A. Knauf and L. C. Lorenz, “Reheating in a Brane Monodromy Inflation Model,” JHEP 0810, 110 
(2008) arXiv:0808.3936 [hep-th]]. 





11 


[24] R. Blumenhagen and E. Plauschinn, “Towards Universal Axion Inflation and Reheating in String Theory,” Phys. Lett. B 
736, 482 (2014) arXiv: 1404.3542 [hep-th]]. 

[25] R. Allahverdi, R. Brandenberger, F. Y. Cyr-Racine and A. Mazumdar, “Reheating in Inflationary Cosmology: Theory and 
Applications,” arXiv:1001.2600 [hep-th], 

[26] M. A. Amin, M. P. Hertzberg, D. I. Kaiser and J. Karouby, “Nonperturbative Dynamics Of Reheating After Inflation: A 
Review,” arXiv: 1410.3808 [hep-ph], 

[27] L. F. Abbott, E. Farhi and M. B. Wise, “Particle Production In The New Inflationary Cosmology,” Phys. Lett. B 117, 29 
(1982); 

[28] A. D. Dolgov and A. D. Linde, “Baryon Asymmetry In Inflationary Universe,” Phys. Lett. B 116, 329 (1982); 

[29] A. J. Albrecht, P. J. Steinhardt, M. S. Turner and F. Wilczek, “Reheating An Inflationary Universe,” Phys. Rev. Lett. 48, 
1437 (1982). 

[30] J. H. Traschen and R. H. Brandenberger, “Particle Production During Out-Of-Equilibrium Phase Transitions,” Phys. Rev. 
D 42, 2491 (1990). 

[31] A. D. Dolgov and D. P. Kirilova, “On Particle Creation By A Time Dependent Scalar Field,” Sov. J. Nucl. Phys. 51, 172 
(1990) [Yad. Fiz. 51, 273 (1990)]. 

[32] L. Kofman, A. D. Linde and A. A. Starobinsky, “Reheating after inflation,” Phys. Rev. Lett. 73, 3195 (1994) |arXiv:hep- 
th/9405187 . 

[331 Y. Shtanov, J. H. Traschen and R. H. Brandenberger, “Universe reheating after inflation,” Phys. Rev. D 51, 5438 (1995) 
arXiv:hep-ph/9407247 . 

[34] L. Kofman, A. D. Linde and A. A. Starobinsky, “Towards the theory of reheating after inflation,” Phys. Rev. D 56, 3258 
(1997) arXiv:hep-ph/9704452 . 

[35] N. McLachlan, Theory and Applications of Mathieu Functions (Oxford Univ. Press, Clarendon, 1947). 

[36] R. Carmona and J. Lacroix, Spectral theory of random Schrodinger operators (Birkhauser, Boston, 1990) 

[37] J. Martin, C. Ringeval and V. Vennin, “Observing Inflationary Reheating,” Phys. Rev. Lett. 114, no. 8, 081303 (2015) 
doi: 10.1103/PhysRevLett. 114.081303 arXiv: 1410.7958 [astro-ph.CO]]; 

[38] J. B. Munoz and M. Kamionkowski, “Equation-of-State Parameter for Reheating,” Phys. Rev. D 91, no. 4, 043521 (2015) 
doi:10.1103/PhysRevD.91.043521 arXiv:1412.0656 [astro-ph.CO]]. 

[39] D. Baumann and L. McAllister, “Inflation and String Theory,” arXiv: 1404.2601 [hep-th]. 

[40] R. Flauger, L. McAllister, E. Pajer, A. Westphal and G. Xu, “Oscillations in the CMB from Axion Monodromy Inflation,” 
JCAP 1006, 009 (2010) arXiv:0907.2916 [hep-th]]; 

[41] H. Peiris, R. Easther and R. Flauger, “Constraining Monodromy Inflation,” JCAP 1309, 018 (2013) arXiv: 1303.2616 
[astro-ph.CO]]; 

[42] R. Easther and R. Flauger, “Planck Constraints on Monodromy Inflation,” JCAP 1402, 037 (2014) arXiv: 1308.3736 
[astro-ph.CO]]; 

[43] T. Kobayashi, O. Seto and Y. Yamaguchi, “Axion monodromy inflation with sinusoidal corrections,” PTEP 2014, no. 10, 
103C01 (2014) arXiv: 1404.5518 [hep-ph]]; 

[44] T. Higaki, T. Kobayashi, O. Seto and Y. Yamaguchi, “Axion monodromy inflation with multi-natural modulations,” 
JCAP10 (2014) 025 arXiv: 1405.0775 [hep-ph]]; 

[45] R. Flauger, L. McAllister, E. Silverstein and A. Westphal, “Drifting Oscillations in Axion Monodromy,” arXiv: 1412.1814 
[hep-th]; 

[46] Y. F. Cai, F. Chen, E. G. M. Ferreira and J. Quintin, “A new model of axion monodromy inflation and its cosmological 
implications,” arXiv:1412.4298 [hep-th]. 

[47] M. A. Amin, R. Easther, H. Finkel, R. Flauger and M. P. Hertzberg, Phys. Rev. Lett. 108, 241302 (2012) arXiv: 1106.3335 
[astro-ph.CO]]. 

[48] S. Y. Zhou, E. J. Copeland, R. Easther, H. Finkel, Z. G. Mou and P. M. Saflin, JHEP 1310, 026 (2013) arXiv: 1304.6094 
[astro-ph.CO]]. 

[49] R. H. Brandenberger, “Lectures on the theory of cosmological perturbations,” Lect. Notes Phys. 646, 127 (2004) hep- 
th/0306071. 

[50] V. Zanchin, A. Maia, Jr., W. Craig and R. H. Brandenberger, “Reheating in the presence of noise,” Phys. Rev. D 57, 4651 
(1998) hep-ph/9709273 ; 

[51] V. Zanchin, A. Maia, Jr., W. Craig and R. H. Brandenberger, “Reheating in the presence of inhomogeneous noise,” Phys. 
Rev. D 60, 023505 (1999) hep-ph/9901207 . 

[52] A. D. Dolgov, A. V. Popov and A. S. Rudenko, “Shape of the inflaton potential and the efficiency of the universe heating,” 
arXiv:1412.0112 [astro-ph.CO]. 

[53] B. A. Bassett, C. Gordon, R. Maartens and D. I. Kaiser, “Restoring the sting to metric preheating,” Phys. Rev. D 61, 
061302 (2000) hep-ph/9909482 ; 

[54] B. A. Bassett and F. Viniegra, “Massless metric preheating,” Phys. Rev. D 62, 043507 (2000) hep-ph/9909353 ; 

[55] F. Finelli and R. H. Brandenberger, “Parametric amplification of metric fluctuations during reheating in two field models,” 
Phys. Rev. D 62, 083502 (2000) hep-ph/0003172 ; 

[56] H. B. Moghaddam, R. H. Brandenberger, Y. F. Cai and E. G. M. Ferreira, “Parametric Resonance of Entropy Perturbations 
in Massless Preheating,” arXiv: 1409.1784 [astro-ph.CO]. 

[57] G. N. Felder and I. Tkachev, “LATTICEEASY: A Program for lattice simulations of scalar fields in an expanding universe,” 
Comput. Phys. Commun. 178, 929 (2008) |hep-ph/0011159 ; 

[58] G. N. Felder, “CLUSTEREASY: A program for lattice simulations of scalar fields in an expanding universe on parallel 




12 


computing clusters,” Comput. Pliys. Commun. 179, 604 (2008) |arXiv:0712.0813 [hep-ph]]; 

[59] A. V. Frolov, “DEFROST: A New Code for Simulating Preheating after Inflation,” JCAP 0811, 009 (2008) 
arXiv:0809.4904 [hep-ph]]; 

[60] J. Sainio, “CUDAEASY - a GPU Accelerated Cosmological Lattice Program,” Comput. Phys. Commun. 181, 906 (2010) 
arXiv:0911.5692 [astro-ph.IM]]; 

[61] R. Easther, H. Finkel and N. Roth, “PSpectRe: A Pseudo-Spectral Code for (P)reheating,” JCAP 1010, 025 (2010) 
arXiv: 1005.1921 [astro-ph. CO] ]; 

[62] Z. Huang, “The Art of Lattice and Gravity Waves from Preheating,” Phys. Rev. D 83, 123509 (2011) arXiv: 1102.0227 
[astro-ph.CO]]; 

[63] J. Sainio, “PyCOOL - a Cosmological Object-Oriented Lattice code written in Python,” JCAP 1204, 038 (2012) 
arXiv: 1201.5029 [astro-ph.IM]]. 



